Realizability of point processes 
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Abstract 

There are various situations in which it is natural to ask whether 
a given collection of k functions, Pj(ri, . . . ,rj), j = l,...,k, defined 
on a set X, are the first k correlation functions of a point process 
on X. Here we describe some necessary and sufficient conditions on 
the pj's for this to be true. Our primary examples are X = M. d , 
X = and X an arbitrary finite set. In particular, we extend a result 
by Ambartzumian and Sukiasian showing realizability at sufficiently 
small densities /?i(r). Typically if any realizing process exists there 
will be many (even an uncountable number); in this case we prove, 
when X is a finite set, the existence of a realizing Gibbs measure with 
k body potentials which maximizes the entropy among all realizing 
measures. We also investigate in detail a simple example in which 
a uniform density p and translation invariant p2 are specified on Z; 
there is a gap between our best upper bound on possible values of p 
and the largest p for which realizability can be established. 



1 Introduction 

A point process in a set X is a random collection of points in X, whose 
distribution is described by a probability measure /i on the set of all possible 
point collections. Here we will take X to be M. d , a lattice such as Z d , the 
torus T d , a lattice on the torus, or an open subset of any of these. We 
always assume that any bounded subset of X contains only finitely many 
points of the collection (this is of course automatically true if X is a lattice); 
the collection of points is then necessarily countable and we will write it as 
{xi, x 2 , . . .}, with the understanding that the Xj are all distinct. 

Well known examples of point processes are Gibbs measures for equilib- 
rium systems of statistical mechanics. The points of the process are then 
interpreted as the positions of particles; because the particle configuration 
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is identified with a subset of physical space, the models satisfy an exclusion 
principle: no two particles can occupy the same position. Point processes 
are also used to model phenomena other than those of statistical mechanics, 
such as trains of neural spikes pO, [2] , departure times from queues [3] , and 
positions of stars in galaxies [I] . 

In many of these cases the quantities of primary interest, partly because 
they are the ones most accessible to experiment, are the low order correla- 
tions, such as the one particle density pi(ri) and the pair density p 2 (ri?i"2)- 
These may be defined in terms of expectations (averages), with respect to 
the measure p, of products of the (random) empirical Geld rj describing the 
process. For continuum systems, i.e., when X is M. d or an open subset of M. d , 
r] is defined by 

»7(r) = (1-1) 

i 

where 5 is the Dirac delta function and the Xj's are as above the (random) 
positions of the points of the process. Then, with (•) denoting expectation 
with respect to /i, 

pi(n) = Mn)), (1.2) 

P2(ri,r 2 ) = (77(^)77^2)) -/>i(ri)<J(ri-r 2 ), (1.3) 
and so on; in general, 

p n ( ri ,...,r n ) = / fl S ( T k-*i k )Y (1-4) 

Equation (jl.4p defines p n as a measure, and we will always assume that this 
measure assigns finite mass to any bounded set in M. dn (which means that 
the number of particles in a bounded set in R d is a random variable with 
finite moments up to order n). In many cases this measure is absolutely 
continuous with respect to Lebesgue measure and we identify p n ( r i> • • • > r n) 
with its density, i.e., p n ( r i> • • • > r n)^ r i • • • dr n is the probability of finding a 
particle in the infinitesimal box dri at for i = 1, . . . ,n. We will assume that, 
whenever possible, p n is extended by continuity to be defined at points where 
two of its arguments coincide. The p n , here and in the lattice case discussed 
in the next paragraph, are often referred to as the n-particle distribution 
functions or correlation functions. 
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When X is discrete, a finite set or a lattice, we will also use the notation 
(11.11) and the definitions (ll.2l) - (ll.4l) . interpreting 8(r — Xj) as the Kronecker 
delta <5r iX i, so that rj(r) has value or 1, depending on whether the site 
r G X is empty or occupied. Note that if n > 2 then p n (i"i, • • • , r n ), G X, 
vanishes whenever = rj for some i ^ j, and that for distinct sites ri, . . . , r n , 
p n ( v i, ■ ■ ■ i r n) is the probability of having particles at these sites. 

In this paper we shall study the following realizability problem: given 
functions pi(ri), p 2 (ri,r 2 ), • • •, Pfc( r i, ■ ■ ■ > r k), nonnegative and symmetric, 
does there exist a point process for which these are the correlations? Since 
only a finite number of correlations are prescribed, the problem may be 
regarded as an infinite dimensional version of the standard truncated moment 
problem [5]. The full realizability problem, in which all the correlations pj, 
j = 1, 2, . . . are given, was studied by A. Lenard [6]. 

Realizability, and the related question of fully describing the realizing 
process, are long standing problems in the classical theory of fluids [7], [8], 
and [9] , recently revived by Torquato, Stillinger, et al. [Till [TT] . This problem 
has been also investigated in the context of the quantum system see e.g. [8], 
|12j . and [13]. An important ingredient in that theory is the introduction of 
various approximation schemes for computing p 2 (ri,r 2 ), such as the Percus- 
Yevick and hyper-netted chain approximations |9J. It is then of primary 
interest to determine whether or not the resulting functions (pi,p 2 ) in fact 
correspond to any point process, that is, are in some sense internally consis- 
tent. If so they can provide rigorous bounds for the entropy of the system 
under consideration. A novel application of the realizability problem to the 
determination of the maximal density of sphere packing in high dimensions 
is discussed in [l4~|. 

Applications of the problem of describing a point process from its low 
order correlations occur in other contexts, for example, in the study of neural 
spikes [U [2] . In this and other physical situations it is natural to consider 
a closely related problem in which the pj, j = 1, 2, . . . , k, are specified only 
on part of the domain X; for example, on the lattice we might only specify 
the nearest neighbor correlations. See [IS] for a similar problem in error 
correcting codes. We will not consider this case further here, except for some 
comments at the end of section 6. 

An important special case is that in which X is M. d , Z d , or a periodic 
version of one of these (a torus), and the point process is translation invariant. 
The specified correlation functions will then also be translation invariant and 



3 



may be written in the form 

Pi(ri) = p, (1.5) 
Pj(r u . . . , Tj) = p j g j (r 2 -r 1 ,...,r j -r 1 ), j = 2,...,k. (1.6) 

As we often work with k = 2, we write g(r) = g 2 (r). We will often state 
our arguments and results in the translation invariant case, but these may 
frequently be extended to the more general situation; when we do not impose 
translation invariance we will use a notation similar to (11.61) : 

3 

Pj(r u . . . ,tj) = JJpi(ri) Gj(r u r 2 , . . . ,r 3 -), j = 2,...,k. (1.7) 

i=i 

We now make some general remarks in order to put the realizability 
problem in context. First, we observe that if the correlations (I1.5D — (jl.6j) 
can be realized for some density p, then they can also be realized, for the 
same functions g 2 ,...,gk, for any p' with < p' < p [16J. To see this, 
note that if r] (r) = S^( r ~~ x «) * s the empirical field with density p, then 
v( r ) = Yl Qi^( r — x «)) where the Qi are independent, identically distributed 
Bernoulli random variables with expectation p 1 j p, is a field with density 
p' having the same value for all gj's. In other words, the new measure is 
constructed by independently choosing to delete or retain each point in a 
configuration, keeping a point with probability p'/p. (We will refer to this 
construction as thinning.) In this light it is thus natural to pose the real- 
izability problem in the following form: given the gj, j = 2, . . . , k, what is 
the least upper bound p of the densities for which they can be realized? It 
is of course possible in the continuum case to have p = oo; for example, if 
gj = 1 for j < 2 < k then for any density p > a Poisson process realizes 
the correlations. For the lattice systems considered here, on the other hand, 
we always have p < 1. 

Lacking a full answer to this question, one may of course ask rather for 
upper and lower bounds on p. A lower bound p > po may be obtained by 
the construction of a process at density p ; we discuss such constructions in 
sections |3HS1 and show that p > under reasonable restrictions on the gj in 
( 11.61) . Upper bounds on p may be obtained from necessary conditions for 
realizability, some of which are described in section [2J 

Beyond the question of realizability one may ask about the number or 
more generally about the types of measures which give rise to a specified set 
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of correlations pj, j = 1, . . . , k. A natural question in the theory of fluids, 
for example, is whether any of these measures are Gibbsian for interactions 
in a particular class; for example, given p\ and p2, is there a Gibbs measure 
realizing these correlations which involves only pair interactions? This ques- 
tion is considered in section [6], where we discuss the nature of the realizing 
measure /i which maximizes the Gibbs-Shannon entropy. 

The outline of the rest of this paper is as follows. In section |2J we discuss 
some necessary conditions for realizability. (In a separate paper [17J we will 
present general necessary and sufficient conditions.) Sections [3JHJ and[5]cover 
proofs of realizability: in section [3] we prove a theorem, a generalization of 
one proven by R. V. Ambartzumian and H. S. Sukiasian [IB], showing that if 
g — 1 is absolutely integrable and g satisfies a certain stability condition then 
the pair (p, g) is realizable for sufficiently small p, with explicitly given higher 
correlations. In section H] we show that the construction can be extended to 
the case in which the third correlation function g% is also specified, showing 
in particular that the realization determined by (p, g) alone is not unique; 
in fact since $3 can take an uncountable number of values there are at low 
values of p an uncountable number of measures realizing (p,g). We note 
also possible extensions to higher order gj. In section [5] we give a variant 
construction for lattice systems, based on the Lee- Yang theorem [T8J [19] . In 
section [6J we show that a problem with specified {p,g), on a finite set, e.g., 
a periodic lattice, may be realized by a Gibbs measure with just one- and 
two-particle potentials whenever p < p. We make some concluding remarks 
in section [TJ and in the appendices discuss in some detail an illustrative 
one-dimensional example and give some technical proofs. 

2 Necessary conditions for realizability 

Clearly, from (11.41) . realizability requires that 




j = 1, . . . , k. 




We also know that the covariance matrix of the field i](r), 




(7 ? (ri)77(r 2 )) - (77(ri)(77(r 2 )) 

P2(ri, r 2 ) + pi(ri)5(ri - r 2 ) - p 1 (r 1 )p 1 (r 2 ) 



(2.2) 
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must be positive semi-definite, i.e., for all functions (p with bounded support, 
j J ip{r 1 )<p{r 2 )S{r 1 ,r 2 )dr 1 dr 2 >0. (2.3) 
If we take ip(r) = lAe ikr for A a bounded subset of X then (12.31) becomes 
/ pi(r 1 )dr 1 + / / e <k ^ ri - ra ^(r 1 ,r 2 )-pi(n)p 1 (r 2 )]dridr 2 >0; (2.4) 

J A J A J A 

conversely, if (12 .4p holds for all k G M d and all A, then (12. 3p holds for all <p. 

In the translation invariant case (12.41) is equivalent to the nonnegativity 
of the infinite volume structure function S(k): 

S(k)=p + p 2 [ e ikr [g{r) - l)dr > 0. (2.5) 

Here we assume J Rd \g(r) — 1| dr < oo; otherwise (12. 5p holds in the sense of 
generalized functions, cf. |20j. There are corresponding conditions on the 
torus T d , the lattice Z d , and the periodic lattice. If equality holds in (12.51) 
for some k then clearly p is maximal: p = p. 

We note also a necessary condition due to Yamada [2T|: if denotes 
the number of of particles in a region Acl, and if 9 is the fractional part 
of the mean of N A , so that (N A ) = k + 8 with k — 0, 1, . . . and < 9 < 1, 
then the variance V A of A^a, 

v a = I I S(v 1 ,r 2 )dr 1 dr 2 

J A J A 

= / pi{r 1 )dr 1 + / / [p 2 (ri,r 2 ) - pi(ri)pi(r 2 )] dridr 2 , (2.6) 
J a JaJa 

must satisfy 

V A > 9(1-9), (2.7) 

because A^ is an integer: 

V A = ((N A -k) 2 )-9 2 = 0(i-0) + (((N A -k)(N A -k-l))) > 9(1-9). (2.8) 

The above necessary conditions all follow from the more general condi- 
tions that we prove in [17J. In summary these say that, given any functions 
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/2( r i) r 2), /i( r ) an d constant f such that, for any n points r 1; . . . , r n in X, 
J2i^j M r u *j) + Ei /i( r i) + /o > 0, we must have 

J J p 2 (r 1 ,r 2 )/ 2 (r 1 ,r 2 )dr 1 dr 2 + J Pl {x)f x {r)dx + f > 0, (2.9) 

AxA 

for all A C M. d . We prove in fact that in the case k = 2, i.e. for the case that 
only pi and p 2 ar e given, f)2.9p is also a sufficient condition for realizability 
under some mild assumptions on the point process. 

We remark that in the case k — 2 all restrictions on p and g beyond those 
arising from nonnegativity of p and of the covariance matrix S of (I2.2p are 
due to the fact that we want r)(r) to be a point process, since we can always 
find a Gaussian process realizing any pi, p 2 with S > [22]. 

We also note that for g(r) < 1 one has 

S(k) > S(0) = lim -±-V A . (2.10) 

A^oo | TV I 

Hence equality in (12.51) implies that the variance V\ is growing slower than the 
volume. Processes with this property are called superhomogeneous and are 
of independent interest, see [231 QUI [21] • As noted above, superhomogeneity 
can hold, for a given g(r), only at the maximal density p = p. 



3 Realizability for small density 

In this section we show the realizability of a given translation invariant p and 
g(r), r G M. d , for sufficiently small p. Our arguments extend immediately to 
the lattice Z d and to the torus, as well as to non-translation invariant pi(r), 
P2( r i? r 2)- Our results are an extension of those given by R. V. Ambartzumian 
and H. S. Sukiasian [IB] , and are based closely on the key idea of that paper: 
given p and g(r) satisfying suitable conditions, one proves the existence of 
a translation invariant process for which the correlation functions p n , n = 
1, 2, 3, . . ., are given by 

p n (rx,...,r n ) = p n Yl a^i-Tj), (3.1) 

l<i<j<n 

and which therefore solves the realization problem for p and g. The ansatz 
f l3.ll) for the dependence on p and g of the higher order correlations (which 
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determines the point process uniquely) is arguably the simplest one possi- 
ble. It corresponds, for n = 3, to the well-known Kirkwood superposition 
approximation in the theory of equilibrium fluids Here, however, we are 
not talking of an approximation to a particular given point process but rather 
of constructing a realizing process or measure p, whose correlations have the 
form (jgZEp . 

To find a point process corresponding to (13.11) . Ambartzumian and Suki- 
asian used the inclusion-exclusion principle which, for any point process, 
relates the correlation functions p n to the probability densities p A (ri, ■ ■ ■ , r n) 
for finding exactly n particles, with positions ri, . . . , r n , in a region A C M d : 

p A (r 1} ...,r n ) 

00 (— l) fc f 

= 2^ — r; / Pn+fc(ri, . ..,r n ,xi,...,Xfc)dxi...dxfc. (3.2) 

Inserting the ansatz (13.11) in (13.21) . one finds an expression for the proposed 
densities: 

p A (r u ...,r n ) = p n H gfa-Ti) 



l<i<j<n 



00 

X 



00 (-n) k f 

i n #( x i~ x ;) n ^( x i - ^) rfx i • • - rfx fc- (3.3) 

fc=0 ' J Ak l<i<j<k i=l,...,n 

j=l,...,k 

It remains to verify that (13.31) in fact defines the probability densities of a 
point process. 

First, note that the quantities Pn( r i, ■■■, r n) are well defined by (13. 3p for 
any value of p whenever there is, for every region A, a constant Ma such that 
(13.11) satisfies 

\p n (r 1 ,...,r n )\<Ml (3.4) 

for ri, . . . , r n in A. The condition (13 .4p is easily verified for many g (see, 
e.g., Theorem 13.11 below). The remaining problem is to prove that the p A 
are all nonnegative. If this is done, then in each region A the collection 
n = 1,2,..., determines a measure p A defining a point process; if A C A' 
then /i A and p A are compatible and by general arguments (Kolmogorov's 
projective limit theorem) there exists an infinite volume realizing measure p. 

Ambartzumian and Sukiasian considered only the case g(r) < 1, r e 
M. d . For this case they constructed a cluster expansion of the Penrose- Ruelle 



8 



type and obtained inequalities of the Groeneveld-Lieb-Penrose type to show 
nonnegativity of each term in a reorganized expansion of the p^. In order 
to extend their result to g's which can be bigger than one, when the cluster 
expansion is no longer positive term by term, we need to use a different 
approach. Recall the definition of the standard grand canonical partition 
function, in the region A, of a particle system with fugacity z, one-particle 
potential V (1) (y), pair potential V(r), and inverse temperature (3 = 1 



-A 2 



^ <-l<i<k Ki<j<k -I ' 



Then (13.31) takes the form 



i>" n tin-*) 

l<i<j<n 



E A (-p,V^\V), 



(3.6) 



with V(t) = - log(0(r)) and V^(y) (= V® (y; r x , . . . , r n )) = J2tt V(y-Ti). 
Note that in (13.61) the one-particle potential V^ l \ and hence also the parti- 
tion function E,\(z,V^ 1 \V), is different for each p^(ri, . . . , r„), depending 
explicitly on n and on the particle positions r 1; ...,r n . The condition (13.41) 
implies that Ha(z, V) is an entire function of z. 

Suppose now that log Ha (z, V^- 1 ', V) is analytic in z in some domain VL 
containing the origin; then Ha(z, V) can not vanish in f2. In particular, 
if (a, b) is the largest interval on the real axis which contains the origin and is 
contained in fi, then H A (z, V^, V) > for a < z < b, since H A (0, V {1 \ V) = 
1. We will apply this observation by finding a domain Q — a disk centered 
at the origin, of radius R — such that for all A, all n, and all ri,...,r n , 
log Ha (z, V) is analytic in f2; then since g(r) > 0, all p£(ri, . . . , r n ) will 
be nonnegative for < p < R. 

These considerations lead to our main result. We write 



C(g) 



/ \g(x)-l\dx. 



(3.7) 



Theorem 3.1 Let g be a non-negative even function on M. d , and suppose 
that (i) C{g) < oo and (ii) there exists a constant b, 1 < b < oo, such that 
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for all n > 1, 

n 

Y[g(x -^)<b (3.8) 
i=i 

whenever x , x 1; . . . , x n satisfy Yli<j #( x i ~~ x j) > 0- Then $-4\) ^ s satisfied, 
and (p, g) is realizable, for all p satisfying 

0<p<(e6C(^))- 1 . (3.9) 

For completeness we state the analogous result on the lattice. 

Theorem 3.2 Let g be an even non-negative function on Z d , and suppose 
that C(g) := J2 x ez d \d( x ) — 1| < oo. Let b be a constant, with 1 < b < oo, 
such that YYi=i9( x i) — b whenever x\,...,x n satisfy Yii<j 9( x i ~ x j) > 0- 
Then (p,g) is realizable for all p satisfying < p < (ebC(g))~ 1 . 

Remark 3.3 (a) The fact that (13.41) holds under hypothesis (ii) of the the- 
orem, with constant M\ = pb 1 ^ 2 independent of A, is immediate. In the 
language of statistical mechanics, this says that the interaction V is stable. 

(b) If g < 1 then hypothesis (ii) holds with 6 = 1, and we recover the result 
of [IB]. 

(c) Hypothesis (ii) also holds if there exists (I) a D > such that g(r) = 
when |r| < D, and (II) a nonnegative decreasing function ip on [D, oo), 
satisfying J °° ^V(i) dt < oo, such that (g(r) - 1) < ^(|r|) [19J. In the 
language of statistical mechanics, (I) says that V^(r) has a hard core; (II) 
says that V(r) is lower regular [25j. In this case one easily obtains an explicit 
possible value for the constant 6. 

(d) Note that, despite the use of results from equilibrium systems, there is no 
reason to expect the realizing measure p giving rise to the p n of (13.1 p to be a 
Gibbs measure with pair potential (unless g(r) = 1, in which case V(r) = 
and p corresponds to a Poisson process). 

Proof of Theorem I3.lt Denote by the m th correlation function for a 
grand canonical ensemble in A with pair potential V, one-particle potential 
inverse temperature (3 — 1, and activity z; as for the partition function 
( 13.51) for this system, these correlation functions depend through on 
r 1; . . . , r n . By the above remarks it suffices to establish analyticity of in 
a disk \z\ < R with R = (ebC{g))~ l , because from 

\og(E A (z,V^,V)) = - I ^(x)dx (3.10) 

az z J A 
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it follows that logE^z, V) is also analytic in this disk. 

To establish the analyticity of k A we in fact show analyticity of all k A L ; we 
proceed as in the classical proof, following in particular section 4.2 of [19]. In 
this proof the Kirkwood-Salsburg equations for the correlation functions are 
written in an appropriately chosen Banach space in the form k A = zip + z%k A 
for some operator % and fixed vector ip. One shows that \\z%\\ < 1 when 
\z\ < R, so that I — z% is then invertible via a power series in z, and a unique 
solution, analytic in z, exists. The primary change in the proof required in 
our case is that one must introduce a dependence of the operator % on the 
sites ri, . . . , r n . We leave the details to appendix [Bj which is probably best 
read with [T5] in hand. ■ 

We next state a generalization of Theorem l3.1l to non-translation invariant 
systems; the proof is omitted. Let X C M d be open, and recall the notation 
p 2 (x,y) = Pl (x) Pl (y)G 2 (x,y) of flU}. 

Theorem 3.4 Let pi and G 2 be non-negative functions on X and X x X , 
respectively, with G 2 symmetric, and suppose that there exists a constant b, 
with 1 < b < oo, such that for all n > 1, n™ =1 G 2 (x , x.j) < b whenever 
x ,xi,...,x n G X satisfy Yl i<:j p 2 (xi - Xj) > 0. Then the pair (pi,p 2 ) is 
realizable if 

eb sup ( / |G 2 (x,y) - l|p x (y)dy ] < 1. (3.11) 

xGX \Jx J 

3.1 Decay of correlations 

We are interested in the decay of the truncated correlation functions Uk for 
the realizing measure specified by (j3.ip . defined recursively by [19J 

n k 

p n (r 1 ,...,r n ) = ^2 Yl Yl u \ii\(( T i)ieij)> ( 3 - 12 ) 

k=l {h,...,I h }eV h (n) j=l 

where t S > } t {n) denotes the set of all partitions of {1, . . . , n} in k disjoint sets. 
We consider only the case in which X = M. d and p\ and p 2 are translation 
invariant; then U\{j\) = p and M 2 (r 1 ,r 2 ) = p 2 [g(r 1 — r 2 ) — 1]. For the 
correlation functions (13. ip the corresponding truncated correlation functions 
have the form 

Mn (r 1 ,...,r n )=p™ Yl II {9(*i-*j)-l)> ( 3 - 13 ) 

GGS c (n) {i,j}£G 
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with S c (n) the set of all connected subgraphs of the complete graph with 
vertex set {1,2,..., n}. 

Let 7(n) denote the set of all undirected trees on {1, . . . , n}. Then from 
(I3.13P and an estimate of Penrose [26J, 

Wr : ,..,r n )| <p n b n - 2 II k^-r.O-ll, (3.14) 

TeT(n) {i,j}eT 

where b is defined as in Theorem 13.11 Using |T(n)| = n n ~ 2 we then obtain 
the L 1 decay property 

/ \u n+1 (r , n,...,p B )| d n v < p n+1 ((n + l)b) n ~ l C (g) n . (3.15) 

Compare Theorem 4.4.8 of [19|. 

One may also establish a pointwise decay bound: if \g(r) — 1| decays 
polynomially or exponentially, then u n (r 1 , . . . , r n ) also decays polynomially 
or exponentially, respectively, with maxi<j <:7 -< n | — rj\. For example, if 
|<?( r ) _ 1| < -Die~ D2 ' r ' for some Di, D 2 > and all r, then 

K(ri r n )| < (^r^rVe-^, (3.16) 

where L is the minimal length of a tree connecting all points ri, . . . , r„ and 
the length of a tree T is ^2nj\ eT — r j | . 

This decay implies that the realizing measure is mixing and therefore 
ergo die. 

4 Triplet correlation function 

We now consider briefly the application of the ideas of section [3] to the real- 
ization problem under the specification of pi, . . . , p/~ for k > 2. For simplicity 
we discuss only the case k = 3 and restrict our considerations to translation 
invariant correlations in M. d which have densities with respect to Lebesgue 
measure. Other cases can be treated analogously. We adopt the notation 

P3(x,y,z) = p 3 g 3 (y - x,z - x) 

= p 3 #(y-xMz-xMz-y)(? 3 (y-x,z-x); (4.1) 

the first equation here is just (jl.6p . and the second is justified by the fact 
that, from the definition (11.41) . p 2 ( x ; y) cannot vanish on a set S of positive 
(Lebesgue) measure unless p 3 {x, y, z) vanishes, for almost all z, if (x, y) 6 S. 
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In analogy with (13.1 j) we make now the ansatz 



p n (r 1} . . . , r n ) := p n JJ #(r; - r^) JJ g 3 (r k - r u tj - r;)(4.2) 

l<«<j<™ l<i<j<k<n 

for the higher correlation functions (n > 4). As before, the probability 
densities p!^ for the point process defined by (14.21) can be written in terms of 
the correlations via the inclusion-exclusion principle (13. 2p and thus in terms 
of a grand-canonical partition function E,\(z, V^ 2 \ V^) for a particle 
system in A with fugacity z, one particle potential V^ l \ non-translation- 
invariant pair potential V^ 2 \ and translation invariant triplet potential V^: 

Pn( r u • • • , r n ) = p n Yl #( r i - r i) II ~ g '^ r i ~ ri ' Tk ~ r ^ 

l<i<j<n l<i<j<k<n 

xE A {-p,V<»,vW,V®), (4.3) 

where 

3 )(x,y,z) := -ln(s 3 (y-x,z-x)), (4.4) 

V^(x,y) : = -\n(g(y-x)) + £ (4.5) 

l<i<n 

V«(x) := ^^ 2 )(x-r,)+ £ ^ 3 )(x, ri ,r,). (4.6) 

l<i<n l<i<j<n 

In order to proceed as in section [3] we have to show that there exists a domain 
for z, independent of ri, . . . , r n , in which In Ha (2, V^ 2 \ V^) is analytic. 
Now, however, we must work with the cluster expansion for both the pair and 
the triplet interactions, as in [27]. We give a result in which the hypotheses 
have been chosen to keep the proof simple and are thus far from optimal. 
We let Vd denote the volume of the sphere in M. d of diameter 1 and write 

C 3 (g 3 ) = sup max{|^(x,y)-l|,|^(x,y)-l| 1/3 }. (4.7) 

x,yGR d 

Proposition 4.1 Let g and g 3 be be given, and assume that (i) g satisfies 
the conditions of Remark \3. 3( c), and C{g) < oo; (ii) there exists a D 3 > 
such that <73(x 2 — x 1; x 3 — Xi) = 1 if |xj — Xj| > D 3 for some i,j,k, and 
C3Q73) < oo. Then {p,g,gs) is realizable whenever 



-l 

< p < 



3 (1 + bC 3 (g 3 )f Ds/D)M (C(g) + v d (D 3 /2) d C 3 (g 3 ))] " , (4.8) 
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where b is defined as in Theorem \3.1\ and b 3 is a constant such that 



IT II ^ 3 ( x * ~~ x °' X J ~ x °) - bs ( 4 ' 9 ) 

j=i i=j+i 

for all n and all x , x 1; . . . , x n with |xj — Xj| > D for < i < j < n. . 

Remark 4.2 (a) Without loss of generality we may assume that D 3 > D, 
since otherwise (73 = 1 and the proposition reduces to Theorem 13.11 
(b) Since there can be at most 

points within a distance D 3 of Xq, all separated from xq and from each other 
by a distance at least D, we may in particular take 

h = (1 + C 3 (S3)) (3D3/D)2d • 

Note that with this choice the upper bound of (14.81) converges as 6*3(^3) \ 
to the upper bound of (13. 9p . 

Proof of Theorem I4.lt The proof is a fairly straightforward extension of 
the proof of Theorem 13 . 11 We sketch some details in appendix [B] ■ 
Finally we would like to point out a consequence of Proposition 14. II which 
illustrates the fact that the same (finite) family of correlation functions may 
be realized by distinct measures, and in fact by mutually singular measures, 
where two measures are called mutually singular if configurations of points 
typical for one of the realizing point process are atypical for the other one, 
i.e., if there exists a set of point configurations A such that A has probability 
1 for one measure and probability for the other. 

Corollary 4.3 Let g be a function fulfilling the conditions of Remark VJ. 3( c). 
Then for any p satisfying the bound Ii3. 9\) of Theorem \3.1\ with strict inequality 
there exist uncountably many distinct and in fact mutually singular realiza- 
tions of (p, g) by point processes. 

Proof: Under the hypotheses one may choose g 3 quite arbitrarily, subject 
only to a condition that C 3 (^ 3 ) be sufficiently small (how small depends 
on -D 3 ), and still have p satisfy the bound (14. 8p . Thus there are certainly 



14 



uncountably many realizations with distinct three point functions. To show 
that these are mutually singular one first establishes, following the procedure 
of subsection 13.11 the decay of the truncated correlation functions for the 
measure constructed in Proposition 14.11 A direct consequence of this decay 
is that the corresponding realizing point process is mixing and therefore 
ergodic. Since any two translation invariant ergodic measures are either 
identical or mutually singular, the result follows. ■ 



5 The Lee- Yang approach 

While for given pj we cannot in general improve the region of realizability 
beyond that described in section [3] and HI there are special situations in which 
more can be said. One class of examples is treated in appendix |A] In this 
section we consider a lattice gas on a countable set X, e.g., X = Z d , with 
G2( x , y) > 1 for x, y 6 X, x ^ y. This enables us to to establish realizability 
using techniques developed for proving the Lee- Yang theorem |18j . 

Theorem 5.1 Let X be a countable set and suppose that (^(x, y) > 1 for 
all x, y G X with x^y and that 

b := sup | [ G 2 (x, y)<oo. (5.1) 

xeX ye X\{x} 

Then (p, G2) is realizable for all 

< p < r 1 . (5.2) 

The condition (I5.2p improves the result of Theorem I3.2[ increasing the 
upper bound on p by a factor of eC(g), with C(g) = X] x ez d \d( x ) ~ M as m 
that theorem. Note that C(g) > 1 since g(0) = 0. 

Proof: We again use the ansatz (13.1 p for the higher correlation functions. 
Let A be a finite subset of X; then (15. ip implies that p n fulfills the bound 
(13.41) with Ma = pb 1 ^ 2 . Hence we can write the probability densities p!^ in 
terms of the correlation functions via (I3.2p . Since the correlation function 
p n+ fc(ri, . . . , r n , xi, . . . , Xfc) vanishes when any of its arguments coincide, one 
can work with the variables £ := {r 1; . . . , r n } and 7 = {x 1; . . . , x^}, where £ 
and 7 vary over all finite subsets of A with 7 fl £ = 0. Then (13. 2p may be 
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written in terms of 7 and £ as = p n (£)H„(£), where 

= E H0 W II G ^y) 1/2 H G ^y) (5-3) 

x e7\{y} y^7 

(compare (I3.5p - fl3.6p ). As before the main problem is to verify that > 0. 
To apply techniques used for the Lee- Yang theorem we write in terms of 
the set a — (A \ £) \ 7 of empty sites rather than in terms of 7. Writing 
A := A \ £ we obtain 



sj(o = (-p) ,Ai n n G2(r ' y) n ^(^y) 172 ) (5-4) 

yeA \ r e£ xeA\{y} 

Ellf-^II^yr 1 ]! ^(^y)- 1/2 n G^y)-^ 

o-cAy e,T \ r 6? xeA\{y} xeA\o- 

Clearly the prefactor here is non- negative in general and is positive for p > 0. 
To prove that the sum is non-negative we rewrite it in the form 

EnU n ^y)> ( 5 - 5 ) 

o-cAyeo- y xeA\o- / 

where 

% := -^n^'.y)" 1 II G a( x .y)" 1/a . ( 5 - 6 ) 

re? xeA\{y}) 

A x , y := G^x.y)- 1 / 2 . (5.7) 

Note that from G2( x , y) > 1 for x 7^ y it follows that — 1 < A x y < 1. Then 
Proposition 5.1.1. of [TjJ] implies that (15. 5p is not zero if \z y \ > 1 for all 
y G A \ f . |zy| can be bounded below by p" 1 rixeA\{y} G ^ y)' 1 > {p^Y 1 ■ 
We have thus shown that H^(£) has no zeros for < p < 1/b. But H^(£) = 1 
for p = 0, so that S^(£) and hence p^(£) is non-negative for all < p < 1/b. 
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6 Gibbsian measures 



In this section we ask whether a specified set of correlation functions pj, 
j = 1, . . . , k, which can be realized by at least one point process, can also 
be realized by a Gibbs measure involving at most /c-particle potentials. Here 
we will first consider this problem for the case in which our system lives on 
a finite set A, e.g., a subset of the lattice. On a finite set every measure is 
Gibbsian in a general sense, so the important restriction is to be Gibbsian for 
a set of potentials involving at most k particles: we will say that a measure 
v on {0, 1} A is k-Gibbsian if it has the form 

1/(17) = z- 1 exp j - ]T ■■■> x ;M x i) • -iM I > C 6 - 1 ) 

y j=i xi^xa^-^Xj-eA J 

where r\ G {0, 1} A , Z is a normalization constant, and —00 < (jfi' < 00. 
As in ( I1.7P we write, for j = 2, . . . , k, 

j 

/ty(xi,...,Xf) = JJpi(xj) G i (xi,...,x J -), xi, . . . ,Xj <G A, (6.2) 

i=l 

and again think in terms of specifying the Gj, j = 2, . . . k, and asking for what 
densities pi(x) the correlations (16.21) may be realized by a fc-Gibbs measure. 
We will prove that this is possible whenever pi(x) satisfies pi(x) < pi(x) for 
all x G A, with equality allowed if pi(x) = 0, for some pi with the property 
that pi and the Gj, j = 2, . . . , k, are realizable. The proof is presented only 
for k = 2, but the result could easily be extended to general k. 

The key ingredient in the argument is the fact that Gibbs measures are 
those which maximize the Gibbs-Shannon entropy of the measure p, 

SO) = -^2n(rj) log nin) (6.3) 
1 

subject to some specified constraints [19J. In particular, if one can use the 
method of Lagrange multipliers to find a measure which maximizes the en- 
tropy, subject to the constraint of a given pi and p2, then the maximizing 
measure will be 2-Gibbsian and the Lagrange multipliers obtained in this 
way will be the desired one body and pair potentials [29]. Here we verify 
that if pi(x) < pi(x) (with equality allowed if pi(x) = 0, as described above) 
then the method of Lagrange multipliers will indeed apply. 
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Theorem 6.1 Suppose that the pair (p 1; G 2 ) is realizable on A. If p\ satisfies 
< Pi(x) < Pi(x) for all x6 A, with pi(x) < pi(x) unless pi(x) = 0, then 
(pi, G2) is realizable by a 2-Gibbsian measure v for some uniquely determined 
potentials 0W(x) ; x G A and (p^(x,y), x, y G A,x 7^ y. Moreover, v 
maximizes the Gibbs-Shannon entropy k6. 3\) over all measures p realizing 
(pi,G 2 ). 

Proof: We begin with a preliminary remark. Suppose that we have verified 
the theorem in the case in which all pi(x) (and hence also all pi(x)) are 
strictly positive. Then the case in which, say, pi(x) = for x G A' C A, is a 
direct corollary: we obtain immediately a 2-Gibbsian measure on {0, 1} A \ A 
realizing the correlations there, and then take </>W(x) = 00 for x G A'. 
Similarly, if G 2 (x, y) = for some pair of sites x, y G A, x 7^ y, with pi(x) 
and pi(y) nonzero, then we set 0^ 2 ^(x, y) = 00, which guarantees that if v 
is given by (16.11) then v{ri) = whenever 7/(x) = i](y) = 1. Thus in the 
remainder of the proof we will assume that pi (x) > for all x G A and prove 
the existence of a realizing measure u, of the form (16. ip with finite potentials 
and with k = 2, on the set of configurations 

C G2 := { V I 77(x)77(y) = if G 2 (x, y) = 0, x + y} . (6.4) 

We now turn to the main body of the proof. For any r] we let |?7| = 
^ xeA ?7(x) be the number of particles in the configuration rj. We first show 
that there exists a measure p* realizing (pi, G2) for which p*(r]) > whenever 
rj G Cg 2 and \r]\ < 2. By hypothesis there exists a measure p realizing 
(pi, G 2 )- We may thin this measure as in section [fl deleting a particle at the 
site x with probability 1 — p 1 (x)/p 1 (x), independently for each site, to obtain 
a measure p* realizing (pi, G2). Now we observe that if i] G Qg 2 an d M < 2 
then p*(^) > 0. For example, if \r)\ = 2 with r/(x) = r](y) = 1 for x^y, 
then G 2 (x, y) > by (16. 4p and hence, since p realizes (pi,G 2 ), pifl) > at 
least for one fj with fj (x) = ^(y) = 1, and there is a positive probability that 
r) will result from fj applying the thinning process. 

Next we construct a measure p realizing (pi,G 2 ) for which /}(//) > for 
all rj G Cg 2 . We first fix e > and for |^| > 2 define p(r/) = p*(^) + e. Now 
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the condition that fx realize {p\,Gq) is that 

E AW = 1, (6-5) 



E 



E 



v(x)Mv) = Pi( x ) ; x G A, (6.6) 
V(*)v(y)fr(v) = Pi(x)pi(y)G 2 (x,y), x, y e A, x ^ y. (6.7) 

Equations (I6.5l) - (l6.7p may be regarded as a system of linear equations for the 
(as yet) undefined fi(f]), \r]\ < 2; note that the number of these unknowns is 
the same as the number of equations. The coefficient matrix in this system is 
(after an appropriate ordering of the rj, \r)\ < 2) upper triangular, with unit 
diagonal; thus these equations can be solved uniquely for the p(r]), \r]\ < 2, 
in terms of these given fi(rj), \r)\ > 2. The resulting fx will differ from //* by a 
perturbation of order e; in particular, since Li*(rf) > for \rj\ < 2, we can by 
choice of e guarantee that also fi{rf) > for \r]\ < 2. But fi{rf) > for > 2 
by construction, so that fx has the desired properties. 

Finally we show that the Gibbsian measure we seek is the measure which 
maximizes S(fx) (see (16.31) ) among all measures realizing (p 1 ,G 2 ). Let v be 
such a maximizer; v is unique by the strict concavity of S. We first observe 
that v must lie in the interior of Cg 2 , i.e., that v{rf) > for all r] e Cg 2 ; 
otherwise define v t = (1 —t)v + tv and note that then 

j t S{u t ) =-^log^) = oo, (6.8) 

so that S(v-t) > S(u) for some t > and v cannot maximize S. Hence v may 
be obtained by the method of Lagrange multipliers, with ( 16.5D — ( IU77T) (written 
for v rather than fx) as constraints. A simple computation shows that the 
v{rf) then have the form (16. ip with k = 2, where the 0^(x) and 0*- 2 ^(x, y) 
are the Lagrange multipliers associated with (16.61) and (16. 7p . respectively (the 
multiplier for (16. 5ft is related to the factor Z). 

To verify uniqueness of the potentials, note that any 2-Gibbsian measure 
realizing (p 1 ,G 2 ) with, say, potentials ip^ 1 ' and tp^, must satisfy the La- 
grange multiplier equations with these potentials as multipliers and is hence 
an extremum of the entropy. From the uniqueness of the extremum and the 
non- degeneracy of the constraint equations ( 16.5l) -( |6T7l) it follows that these 
multipliers are uniquely defined, i.e., that ip^ = 0^ and = (p^. ■ 
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6.1 Infinite volume 



It is natural to ask if there exists an analogue of Theorem 16.11 for an infinite 
lattice, such as Z d , for example when the given correlation functions are 
defined from the beginning on Z d as translation invariant quantities, i.e., 
Pi(x) = p, G 2 (x, y) = g(y — x) as in (11.51) . (11.61) ; on Z d we will call a measure 
2-Gibbsian if it satisfies the DLR equations for an interaction with only one 
and two body potentials. A result in this direction is due to L. Koralov [30J; 
using cluster expansion techniques, he has established the existence of an 
infinite-volume 2-Gibbsian measure in the lattice case for k = 2, p small, and 
g sufficiently close to 1 — specifically, for \d( r ) ~ M — 1- 

An attractive alternative approach would be to apply Theorem 16.11 in 
large boxes A C Z d to obtain potentials ( x ) and 0^ (x, y), x, y G A, 
realizing (p, g) in A, and then to show that under suitable restrictions on p 
and g the A /* 7h d limits of these potentials exist and are summable and 
translation invariant. We do not know how to carry out such a program, 
as we have no control on the behavior of the <f)\ as A changes. In fact we 
do not know if the 2-Gibbsian measures u& realizing (p, g) in A converge 
as A /* Z d to any measure v on Z d . Any sequence of such measures must 
have a (weakly) convergent subsequence, however, by compactness, and the 
limiting measure v will realize the translation invariant (p, g) on Z d . In [31] 
we showed that any such realizing measure which is translation invariant and 
Gibbsian, with summable potentials, is necessarily a 2-Gibbsian measure with 
uniquely determined potentials which maximizes the entropy density among 
all realizing translation invariant measures. 

It follows from the above that if there is in fact a translation invari- 
ant entropy maximizing 2-Gibbsian measure v on Z d , realizing (p,g), with 
summable translation invariant potentials ip^ 1 ' and ^^(x — y) is just 
the chemical potential), then the conditional measure ^(-|^a c ) on A, for a 
specified configuration rp^c on Z d \ A, will be the 2-Gibbsian measure which 
maximizes the entropy for (pi(x\rj\c), G^fx, y\r}A c )), the one an d two particle 
distributions obtained from ^(%|^ac). If furthermore this measure is unique 
for these potentials, then clearly iy(-\r]\c) — > v for every rp^ as A /* Z d . 

To obtain translation invariant measures we may take for the domain 
A in Theorem 16.11 a periodic lattice L, where L = {— L + 1, . . . , L} d with 
periodic boundary conditions, and find a measure z/l realizing (p, pg^) for 
some periodic g^(r) defined for those r satisfying |rj| < L, i — 1, . . . , d, and 
such that gh( r ) —* g( r ) as L oo. In this case any subsequence limit v of 
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v%, will be translation invariant in Z d and realize (p, pg) . Conversely, from 
a translation invariant p we can construct Pl by first projecting p into a 
cubical box A of side 2L to obtain p\, then defining, for r] a configuration 
in A or equivalently L, Pl(??) = (2L) _d ^ xgL /i(r x ^), where r is the shift 
operator on L. This yields a periodic measure pl with density p and with 
g L (r) = g(r) + 0(1/ L) for fixed r. 

The real question then is whether any subsequence limit of the z/jl will be 
a 2-Gibbs measure with summable pair potentials. To answer this requires 
some control of the potentials 0V( X — y)> which we lack at present (see 
question 4 at section [7]). 

7 Concluding remarks 

There are clearly many natural questions left unanswered by our results. We 
list some of them for the case of a specified p and g(r). 

1) Is there any practical way to bridge the gap between the obvious nec- 
essary conditions described in section [2] and the sufficiency conditions given 
in sections [3] and \5tf 

2) When is the measure defined by (13.11) Gibbsian or quasi-Gibbsian? 

3) Can one extend Theorem 16.11 to continuum systems in a finite domain 
A C R d ? We expect this to be true under some reasonable assumptions on 
g(r), e.g., the hard core condition that g(r) = for r < D, D > 0, under 
which there can only be a finite number of particles in A. 

4) What can one say about the existence and nature of an entropy maxi- 
mizing measure on Z d for (p, <?)? In particular are there situations for p < p 
when such a measure is not a translation invariant 2-Gibbsian measure? As 
pointed out at the end of section 0, if the answer to this question is no then 
the v obtained from the periodic of formula (16. ip will be the 2-Gibbs 
entropy maximizing measure on Z d . 

5) What happens to the realizability problem if one does not specify g(r) 
for all r 6 Z d (or IR d ), but only for r in some finite domain, say |r| < Rl 
(This question was mentioned briefly in sectionTJ) As the notation indicates, 
we are still considering translation invariant correlations; we may in addition 
require that the two-point correlations p2(x, y) of the realizing measure all be 
translation invariant and approach p 2 as |x — y| — > oo. This is the problem 
discussed in [15] for r G Z, or just on a ring. It turns out that, at least 
for the case r 6 Z and specified (p, g(l), . . . , g(k)), one can compute, via 
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a finite number of operations, whether the correlations are realizable. On 
the d dimensional periodic lattice L, for general d, the entropy maximizing 
measure in (16.11) will now contain only a finite number of terms, so that 
the transition to an entropy maximizing 2-Gibbs measure on Z d with only 
finite range potentials appears feasible. This will be described in a separate 
publication. 



A Simple examples 

In this appendix we collect some realizability results for certain concretely 
given g. Let us first mention a realizability problem in M. d which has been 
extensively studied by Torquato and Stillinger []JJ [TT]: to determine for 
which densities p there exists a translation invariant point process on M d 
with 

ff (r)={°' "H" 1 ' (A.1) 
V ; \l, if |r| > 1. V ' 

Condition (12.51) implies that (p, g) can be only realized if p < (v ,j2*) -1 , where 
Vd is the volume of the ball with diameter 1 in M. d (v\ = 1, t>2 = vr/4, etc.). In 
the other direction, Theorem 13.11 implies that for general d these correlations 
are indeed realizable if p < e~ 1 v d ~ 1 2~ d . Thus the maximum density p(d) for 
which g is realizable satisfies 

er 1 < 2 d v d p(d) < 1. (A.2) 

Ideas introduced in recent work on the radius of convergence for pure hard- 
core systems [32] can be used to give a stronger version of Theorem 13.11 
improving the lower bound in (I A. 21) beyond e _1 , e.g., for d — 1 to 0.4 and for 
d — 2 to 0.5107 (numerically). In one dimension we can say more: a simple 
construction of [33J shows realizability by a renewal process if p < 1/e (it is 
also shown in [3BJ that 1/e is the maximal density for which a renewal process 
can realize g). A more complicated construction p5TJ, using hidden Markov 
processes, gives realizability for all p < 0.395, so that 0.395 < p(l) < 0.5. 
The gap between these upper and lower bounds remains as a challenge to 
further rigorous analysis. Certain simulation results [TT] have suggested that 
in low dimensions the process may in fact have p(d) = 2~ d v d ~ 1 . However, we 
can show that this is false in one dimension, i.e., that p(l) < 1/2; we discuss 
this further below. 
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This continuum problem is, in dimension d = 1, related to the following 
lattice problem: for what densities p can the second correlation function 
p 2 g^ a \ where 

{0, ifx = 0, 
a, if |x| = 1, (A.3) 
1, if \x\ > 1, 

be realized by a point process on Z? From the remarks in section [1] we know 
that for fixed a the set of realizable densities p is an interval [0, pj with 
< p a < 1. There is a superficial similarity between the continuum problem 
( lA.ip and the lattice problem flA.3[) for a = 0, but there is also a deeper 
relation. For suppose that r] c is a point process in M which realizes flA.lj) at 
density p. Then, for any fceZ, define r] k = N( k)k+ i] = J kl 7] c (x) dx (that is, 
r)k is the number of points of the process r] c lying in the interval (k,k + 1]). 
Then r\ k has value or 1, (rjk) = p, and for j > 0, 

(VkVk+j) = (VoVj) = P 2 dx dyg(y-x) 

Jo Jj 



p 2 [ dx f dy = p 2 /2, if j = l, 

y J p, (a.4) 

p 2 dx dy = p 2 , if j > 2. 



o Jj 

Thus f] solves the lattice problem with the same density p and with a = 1/2, 
from which p(l) < p 1 / 2 . We will see below that pi/ 2 = 1/2, so that the 
relation p(l) < py 2 is consistent with the possibility p(l) = 1/2 discussed 
above. 

We now discuss the lattice problem (IA.3I) in some detail, as an illustration 
of the difficulties to face in the general situation. Clearly pi = 1, since for 
each p G [0, 1] the Bernoulli or product measure v p realizes ( IA.3[) : for p ^ 0, 1 
there are in fact uncountably many mutually inequivalent realizing measures 
[3T] , while for p = 0, 1 the realization is unique. For other values of a, upper 
bounds on p a are provided by (12. 5p and (12. 7p . In particular, (12.51) yields the 
upper bound p a < -Rf(oj), where 

i — — , if 1 > a > 0, 
R F (a) = I 3 ~ l 2a (A.5) 

ifa>l. 



2a - 1 
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0.0 0.2 0.4 0.6 0.8 1.0 

a 

Figure 1: Upper bounds Rf(oz) (solid) and Ry(a) (dashes) for p a . Lower 
bounds rs(a) (dots/dashes) and r^(a) (dots), for < a < 1. The inset plots 
the difference Rp(a) — Ry(a) on a logarithmic scale. 



The Yamada condition (12.71) gives an upper bound p a < Ry(a); a straight- 
forward but somewhat lengthy computation shows that Ry(a) = Rf(c() for 
a = 1/2, for a — (k ± l)/2k, k = 1, 2, . . ., and for a > 1. For other values 
of a, .Ry(a) < i?ir(a), so that certainly the bound Rp is not always sharp. 
These bounds, together with several lower bounds for p a obtained below, 
are plotted in Figure 1; for values of a at which Ry(a) < Rf(oc), Ry{oi) is 
determined numerically. 

Lower bounds on p a come, essentially, from procedures for explicitly real- 
izing the desired process at some value of p. For example, from Theorem 13.11 
for a < 1, and Theorem 15. 1[ for a > 1, we obtain p a > r^(a), where 

if 1 > a > 0, 

r A (a) = { e } ~ za ) (A.6) 

if a > 1. 

Comparison with (IA.5H shows that, as might be expected, the lower bounds 
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from these general construction methods do not approach the upper bound 
very closely. To get better bounds or exact values for p a one must turn 
to more ad hoc methods. In this spirit we next describe two families of 
processes which realize (I A . 3 [) and which provide improvements in the lower 
bounds flA.6j) . We studied other constructions, partially improving some of 
the results below, but these are omitted for conciseness. 

The first construction, valid for a > 1/2, achieves a density p = rs(ct), 
where 

f , if 1 > a > 1/2, 

r s ( a ) = t 1 + V2-2Q ( A . 7 ) 

thus p a > rs(oi). Comparison of (1 A. 71) with (1A.5I) shows that rs(a) = Rf{o) 
for a = 1/2 and a > 1, so that p a = Rp(a) for these values. The measures 
for these processes are superpositions of two measures of period two. To 
construct them, we first choose with equal probability one of two partitions 
of Z, either . . . U {-2, -1} U {0, 1} U {2, 3} U ... or ... U {-1, 0} U {1, 2} U 
{3,4} U . . ., and then assign a configuration to each pair (i, i + 1) of sites 
in the partition independently, taking (r]i,r] i+ i) to have value (1,0), (0,1) 
each with probability p, (0, 0) with probability q, and (1, 1) with probability 
(1— p— q)/2. The optimal choices of parameters which lead to (1A.7I) are q = 0, 
p = V2 - 2a/(l + V2-2a) for 1/2 < a < 1 andp = 0, q = (2a-2)/(2a-l) 
for a > 1. 

The measures constructed above, as superpositions of period-two mea- 
sures, do not have good mixing properties (this defect will be inherited by 
measures with lower values of p obtained via the thinning process described 
in section [T]). For the case a = p = 1/2 this decomposability is inevitable; the 
system is then super homogeneous (see section [2} in fact the variance of the 
number of points on any set of consecutive lattice sites is uniformly bounded 
by 3/4) and the decomposability of any realizing measure then follows from 
a result of Aizenmann, Goldstein, and Lebowitz [35J . 

We now sketch briefly the argument, mentioned in the first paragraph 
of this appendix, that yields p(l) < 1/2. Given a point process r\ c in R 
realizing flA.ip at density p, we define a lattice process r/ n \ n > 1, by rj^ 1 = 



25 



f(k-i)/n ^ c i x ) Then r/ n ) has density p/n and two point function 

if \j\ < n, 

i'lu'lj) = { P 2 /(2n 2 ), if \j\=n, (A.8) 
2 /n 2 , if \j\ > n. 

For n = 1 this is just the construction given earlier (see ( 1A.4j) ). and as there 
we see that if p* is the maximum density for which (1A.8I) can be realized 
then necessarily p(l) < np^; thus p(l) < 1/2 follows from p* < 1/4. To 
verify the latter, one first shows that any realizing measure pi for ( 1 A. 81) with 
n = 1, p = 1/2 is supported on configurations such that for some z with 
z = 1, 2, each pair of sites (2j + i, 2j + i + 1), j G Z, contains exactly one 
particle; this follows from the fact that the random variable ~52i=j Vi (i £ ^? 
fcGN) has minimal variance 1/4 and hence takes only the values k and k + 1. 
(Note that these configurations form the support of the a = 1/2 measure 
constructed above.) Such a pi may be obtained by projecting a realizing 
measure p2 for ( 1A.8I) with n = 2, p = 1/2 via either := 77^ + ^fc+i or 

^fe := ^2fe-i + ^2k 1 an< i if f nen follows that P2 is supported on configurations 
with the property that for some i with 1 < % < 4, all sites 4j + i, j G Z, 
are empty and each triple of sites (4j + i + 1, 4j + % + 2, 4j + % + 3) contains 
exactly one particle. This leads to a decomposition p 2 = Ylt=i ^2,1 an d from 
this we derive a contradiction, since the variance with respect to p2,i of one 

of the random variables Yl^i vfk+i^ ^ = 1, • • • , 4, would be negative for large 
enough m. For details see |34j . 

Our second construction is valid for < a < 1. For this process we 
first distribute particles on Z with a Bernoulli measure such that each site 
is occupied with probability A; then if in this initial configuration a site is 
occupied we delete the particle occupying its left neighbor, if it exists, with 
probability k. With the optimal choices A = 1/(1 + y/1 — a) and k = y/1 — a 
we obtain a realization of (1A.3j) with density 

r B (a) = - ■ \ - , (A.9) 



vT- M 



so that p a > Tgicx). For a > 1/2 this is of no interest, since tb < rs, but for 
0<a<l/2it improves on the bounds flA.6j) . 



The case a = merits special discussion. From (IA.9I) . p > 0.25; in 
this case the point process used to obtain rg is a renewal process and the 
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construction is a lattice version of that given in [33] to establish that p > 
1/e for the continuum problem flA.10 . A construction based on a hidden 
Markov process [31] improves this lower bound to po > 0.265. The upper 
bound po < Rf(&) = -Ry(a) = 1/3 can be improved [31] to p < (326 — 
v / 3115)/822 ~ 0.3287. As in the continuum problem flA~T1) . it remains a 
challenge to diminish the rather large gap between these upper and lower 
bounds. 



B Proofs of Theorems 13.11 and 14.1 



Completion of the proof of Theorem 13.11 We must show that the 
functions are analytic for \z\ < R. As indicated in section [3J the proof 
follows closely the proof of Theorem 4.2.3 of [19], and we content ourselves 
with pointing out a few key steps and the necessary changes. Let E% be the 
Banach space of all sequences ((p m )™ =1 , where <p m : (R d ) m — > C, for which 

IMI^ := sup sup £~ m (<p m (xi, . . . ,x m )| < oo, (B.l) 

m>0 xi,...,x„ 



and let Xm( x i, ■ ■ ■ , x m) be the characteristic function of the set 

s m / m n \ v 

(xi, . . .x m ) g (R d r i x 4 g a, n n ^ - x,) n^x, - ^ > o . 

^ 4— 1 \ o— 4 — 1 I ) 



1=1 \j=i+l j=l 

(B.2) 

Define the operator % on |J^>o by 

(X^) m+1 (x 1; . . . ,x m ,x) = Xm+ i(xi, • • . ,x m ,x)e- v(1) W-^ m =i^ (2) (x-x J ) 

oo _ „ k 

fc=max{0,l— m} J=l 

• v? m +fc(xi, . . . , x m , zi, . . . , Zfe) dzi • • • dz fc .(B.3) 

The factor Xm+i in ( IB. 31) is needed in the estimate ( 1B.4I) below because the 
one-body potential in ( IB. 31) depends on ri, . . . ,r n , and it is primarily 
these aspects which distinguish the proof here from that of [19J. Note that, 
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with this factor, % depends on A, on g, and on ri, . . . , v n . Then 
|(^) m+1 ( Xl ,...,x m ,x)| < |M| e e-^*)-£r=x^(x-x 3 .) 



•Xm+i(xi, . • • ,x m ,x) V i / TT | e - y(2) ( x -^ - 1 e^Z! • • -dz k 

< ||^|| 5 r bexp Uj^ \g{z) - l\dz\ , (B.4) 
so that 

\\X<p\\i < T^exp Uj^ \g(z) - l\dz\ |M|<. (B.5) 

Now as in [1QJ the sequence of correlation functions k^ satisfies the Kirkwood- 
Salsburg equation 

k A = zip + z %k A (B.6) 

where ip m = 5 m i] by an optimal choice of £ we may insure that ||^3C||^ < 1 
whenever \z\ < R. For such z, I — z% is invertible and the equation (1B.6I) 
has a unique solution in one proceeds as in [19]. ■ 
Proof of Theorem 14.11 As in the proof of Theorem 13.11 we denote by 
k^ the m th correlation function for a grand canonical ensemble in A with 
activity z, inverse temperature /3, and potentials (I4.4p - fl4.6p ; we do not 
write explicitly the dependence of these functions on r 1; . . . , r n . The k A again 
satisfy the Kirkwood-Salsburg equation (1B.6[) in the Banach space but 



with a modified operator %. To define % we write Uj := for i = 1, . . . , n, 
u, := Xj_ n for i = n + 1, . . . , n + m, and Uj := Zj_ n _ m for i = n + m + 
1, . . . , n + m + k. Then 

oo „ 

(X^) m+1 (x 1 ,...,X m , X ) = ze -^\u 1 ,-,u n+m ) £ / (R7) 

fc=o ' ^ Ak 

■i^(x|lli, . . . , U n _|_ m _|_fc)<y2 m _|_fc(u m+ i, . . . , U n+m+ fc) d\l n+m+ i ■ ■ ■ <iu n _|_ m+ fc, 
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with 



n+m 

£(x| Ul ,...,u n+m ) := V r(2) (x,u,)+ J2 ^ (3) (x,u,, Uj ) (B.8) 

i=l l<i<j<n+m 

K(*\ 

■ s n ^ e - y(2)(x ' u " +m+i) - 1) 

77C{l,...,fc} ie{l,...,fc}V? 

H (' e -^)(x J u ij u,)_ 1 \ J-J e -^)(x,u ft+ro+i )_ ( B 9 ) 

G {i,j}eG i&7 

Here extends over the set of graphs which have vertex set V\ U V2, where 
V\ = {1, . . . , m + n} and V 2 = {l + m + n \ I <E i]} (& graph being identified as 
a set of edges, i.e., of unordered pairs of vertices), and in which every edge 
has at least one of its vertices lying in V2, and every vertex in V2 has at least 
one edge incident on it. The derivation of flB.7HB.9l) is similar to that of the 
usual Kirkwood-Salsburg equation [19J. 

The operator DC can be bounded as follows. First, the hypotheses of 
the theorem imply that Xn+m+Hi( u , x)e~ E ^ u ^ < b b 3 . Next, since the fac- 
tor (e _y(3> ( x,Ui ' u j) — 1) in (1B.9|) vanishes unless Uj and Uj lie inside the ball 



of radius D 3 around x, the sum over rj is nonzero only for those f] 's such 
that all points (z^)^ are inside this ball. On the other hand, the factor 
Xn+m+Hi( u i x ) implies that K vanishes unless all the Uj are separated by a 
distance at least D, so that we may suppose that there are at most iV (see 
(I4.10p ) of these points inside the ball. We thus have the bound 



JA k 

jr^c( g r\v d (D 3 /2yy f; 



< 

1=0 ^' j=\l/2] 

< (l + b C 3 (g 3 )) M ( C (9) + v d (D 3 /2) d C 3 (~g 3 )) k ■ (B.10) 

Here M = N(N — l)/2, C 3 = sup xygR d |g 3 (x, y) — 1|, \s] is the least integer 
not smaller than s, and we have used the inequality C% < C 3 (g 3 yC 3 (g 3 ) 1 , 
valid for j > 1/2, which follows from (14.71) by considering separately the 
cases C 3 > 1 and C 3 < 1. 
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From ffl~TUj) we have that M < (3D 3 /D) 2d , and from flB~7l) and flRTUj) it 
then follows that the norm of % in the Banach space satisfies 

||3C|| { < C l bh (1 + bC(~g 3 )f D3/D)2d e^ + ^ D ^ dc ^). (B.ll) 

An optimal choice of £ again shows that ||^3C||^ < 1 when \z\ < R3, where 
R3 denotes the right hand side of (14.81) . completing the proof. ■ 
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